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ABSTRACT 

Using high signal-to- noise spectra (S/N 30) taken with the Cosmic Origins Spectrograph (COS) 
of the z em =0.9754 quasar PG1148+549, we report on the physical conditions of Ne VIII + OVI 
absorption line systems at z a b s =0.6838, 0.7015, 0.7248. At these redshifts, absorption lines from 
multiple ionization stages of oxygen (O II, O III, O IV, O VI) are available to constrain the ionization 
and physical conditions of the gas. A single-phase model with only photoionization or collisional 
ionization fails to reproduce the observed column density ratios. Even with extreme modifications 
to the ionizing background, single-phase photoionzation models fail to produce enough Nc VIII. We 
thus conclude that the O VI and Ne VIII are produced via collisional ionization. The mctallicities 
of these systems are determined to be [0/H]> —0.5, higher than the typically assumed value for O 
VI absorbers of [0/H]~ —1.0. Densities of n# ~ 10 -4 cm -3 are determined for the photoionized 
phase, indicating that these systems reside in regions of overdensity A ~ 80. The collisionally ionized 
component of the gas traced by O VI and Ne VIII bears temperatures of T~ 10 5 7 K. We determine 
the redshift density of Nc VIII absorbers cW /dz ~ TL$, similar to the redshift density of O VI 
absorbers of comparable equivalent width to these systems. Based on the redshift density and galaxy 
luminosity function at z~0.7, we estimate the cross section radius of the O VI + Nc VIII absorbers 
to be ~ 70-150 kpc. We find a star forming ~L* galaxy at the redshift of the z Q b s =0.7248 system, at 
an impact parameter of 217 kpc. 

Subject headings: Intergalactic Medium 



1. INTRODUCTION 

The cosmic baryon fraction has now been well con- 
strained by both big bang nucleosynthesis coupled with 
deuterium measurements (O'Meara et al. 2006) as well 
as observations of the cosmic microwave background 
(Spergel et al. 2007) to be f b = 0.17±0.01. At odds with 
the cosmological value, the baryon fractions in nearby 
galaxies as measured by their stellar and gas components 
compared to the dark matter mass are known to be well 
below the cosmic value (McGaugh et al. 2010). One ex- 
planation for these "missing" baryons is that low density 
gas residing in the halos of galaxies and groups is heated 
during infall into the potential well of the structure. This 
shock heating is a universal prediction of hydrodynamical 
simulations of structure formation (Dave et al. 2001; Cen 
& Fang 2006; Tepper-Garcia et al. 2011). Alternatively, 
it is possible that the baryons never cool and fall into the 
dark matter halos as galaxies evolve (Mailer & Bullock 
2004; Anderson & Bregman 2010). In either case, the 
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observations indicate that the majority of the baryons 
in the Universe are not located in the disks of galaxies. 
The most sensitive method to search for such low density 
plasmas currently is via the ultraviolet (UV) absorption 
lines imprinted on the spectrum of a background quasar 
(QSO). 

Simulations and analytical models of galaxy formation 
and evolution require massive outflows to account for 
the observed properties of galaxies such as the stellar 
and ISM mctallicities, as well as the total stellar masses. 
Similarly, the "closed box" model of galactic chemical 
evolution dramatically fails to reproduce the distribu- 
tion of stellar metallicities in the Milky Way and nearby 
galaxies (Worthey et al. 1996; Schmidt 1963). Feedback 
processes from current or recent star formation, cither 
through radiatively driven or supernova driven galactic 
scale winds are thought to play crucial roles in the evo- 
lution of a galaxy. Such outflows are commonly seen in 
nearby star forming galaxies (Heckman et al. 1995; Mar- 
tin et al. 2002; Veilleux et al. 2005), and are thought to 
be responsible for injecting energy and mass into the cir- 
cumgalactic medium (CGM) and beyond into the IGM. 

The O VI AA 1031.9, 1037.6 lines arising from 5+ , 
which has a peak ionization at ^300, 000 K in collisional 
ionization (Gnat & Sternberg 2007), have been used to 
trace hot gas in a range of environments from the local 
ISM to the IGM (Bowcn et. al. 2008). Recent obser- 
vations indicate that the extended halos of star forming 
galaxies are ubiquitously filled with oxygen and nearly 
all have detectable O VI absorption in the halos out to 
~150 kpc. This circumgalactic material appears to con- 
tain a substantial amount of mass, comparable to the 
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Table 1 

Cosmic Origins Spectrograph Observations of PG1148+549 



COS Grating 


Acen 


Observation 


Exposure 


MAST ID" 






Date 




Time 






(A) 






(seconds) 




G130M 


1309 


2009 Dec. 


30 


3192 


LB1022Q3Q 


G130M 


1309 


2009 Dec. 


30 


3192 


LB1022RAQ 


G130M 


1309 


2009 Dec. 


30 


2527 


LB1022PMQ 


G130M 


1327 


2009 Dec. 


26 


2527 


LB1023EBQ 


G130M 


1327 


2009 Dec. 


30 


3192 


LB1022RHQ 


G130M 


1327 


2009 Dec. 


30 


3192 


LB1022RRQ 


G160M 


1600 


2009 Dec. 


26 


3192 


LB1023EOQ 


G160M 


1600 


2009 Dec. 


26 


3192 


LB1023EFQ 


G160M 


1600 


2009 Dec. 


26 


3192 


LB1023EVQ 


G160M 


1623 


2009 Dec. 


25 


2482 


LB lO 24 ACQ 


G160M 


1623 


2009 Dec. 


25 


3192 


LB1024AFQ 


G160M 


1623 


2009 Dec. 


26 


3192 


LB1023F2Q 



of oxygen are well as many other ions are available at 
600-1000 A in the rest frame (Verner ct al. 1994). At z ~ 
0.5, these transitions enter the bandpass of the Cosmic 
Origins Spectrograph (COS) onboard the Hubble Space 
Telescope (HST). 

Aside from O VI, other ions of highly ionized species 
exist with absorption lines in the FUV. The Ne VIII AA 
770,780 doublet enters the COS bandpass at z ~ 0.45. 
The ionization potential of Ne +6 is 207 cV; as such Ne 
VIII is expected to be a more reliable tracer of collision- 
ally ionized gas. Several absorption systems bearing Ne 
VIII have been discovered, indicating highly ionized and 
multiphase gas with temperatures of T=10 K. (Sav- 
age et al. 2005; Narayanan et al. 2009, 2011; Tripp et al. 
2011). 

The excellent sensitivity of COS in the UV gives the 
opportunity to study the IGM and CGM in much more 
detail than has ever been achievable in the past. In this 
work we study in detail absorption systems at z ~ 0.7 
where we have access to a number of ionization states 
of oxygen (O I to O VI) and the Ne VIII AA 770, 780 
doublet towards PG1148+549. The organization of the 
paper is as follows: In § 2, we discuss the COS observa- 
tions and data handling, details of the analysis of each 
absorption system are given in § 3, § 4 describes the ion- 
ization modeling, results are given in § 5, a discussion of 
the galaxies in the field are given in § 6, and § 7 gives a 
summary and discussion. 

2. OBSERVATIONS 

The observations presented here were taken as part of 
program G011741 (PI Tripp), a blind survey of z ~ 1 
UV-bright quasars in search of highly ionized species such 
as O VI, Ne VIII, and Mg X as well as the multiphase ha- 
los of galaxies. Targets for this sample have been selected 
only based on their redshift (z > 1) and flux in the far 
ultraviolet (FUV). PG1 148+549 was previously observed 
with the Faint Object Spectrograph (FOS, R ~ 1300) in 
programs 4952 and 6210, although the resolution and 
S/N are far too low for the purposes of the analysis in 
this work. These COS spectra were acquired on 25-30 
December 2009 and cover a wavelength range of 1150- 
1800 A. Details of the COS observations are given in 
Table 1. 

These data were processed in the same manner as de- 
scribed in Meiring et al. (2011). Multiple focal-plane 
positions were used to reduce the effects of fixed-pattern 
noise when the various exposures were combined. In ad- 
dition, flatfielding was applied. The reduced and coad- 
ded spectra have been binned by 3 pixels, as the raw COS 
data are ovcrsampled with a ~ 6 pixel wide resolution el- 
ement. All subsequent measurements and analysis were 
performed on the binned spectra. The binned spectra 
have a resolution of ~15 km s _1 per resolution element 
and have been normalized using a cubic spline function 
to fit the continuum. The S/N of these COS spectra 
range from 20 — 40 per resolution element. 

3. NOTES ON INDIVIDUAL SYSTEMS 

A number of absorption line systems are observed 
in the spectrum of PG1148+549. Here, we focus on 
three systems at z = 0.6883,0.7015,0.7248 with Ne 
VIII detections as this redshift allows for the multi- 
ple ionization stages of oxygen as well as the Ne VIII 



a Exposure identification code in the Multimission Archive at 
Space Telescope (sec http://archive.stsci.edu/index.html). 



mass of the ISM in the galaxy itself (Tumlinson ct al. 
2011b; Tripp ct al. 2011; Prochaska et al. 2011). 

Hydrodynamical simulations of galaxy evolution show 
that O VI absorbers typically reside in metal enriched 
regions with overdensitics A = p/(p) of 1 — 100 (Ccn & 
Fang 2006; Dave et al. 2001; Tepper-Garcia et al. 2011; 
Oppenheimer & Dave 2009; Oppenheimer et al. 2011; 
Smith et al. 2011). Simulations by Tepper-Garcia et 
al. (2011) indicate that O VI does appear to trace pri- 
marily shock heated gas, with the distribution of tem- 
peratures peaked at ~ 10 5 3 K. However, ^30 percent 
of the O VI systems in that simulation are at lower 
temperatures, and primarily ionized by the local UV 
background. Other simulations indicate that the O VI 
and even Ne VIII absorbers arise primarily in photoion- 
ized gas (Oppenheimer & Dave 2009; Oppenheimer et 
al. 2011). Interestingly, ~ 30% of the observed O VI 
absorbers have characteristics of cool, photoionzed gas 
(Tripp et al. 2008). 

Despite recent advances in the UV spectroscopy of the 
low-z IGM, several questions are still unanswered due in 
part to the lack of diagnostic lines at low redshift. The 
baryonic content is still largely unknown, as well as the 
metallicity and ionization mechanisms of the gas com- 
prising the IGM and CGM. Much of the work on \ow-z 
absorption lines has been accomplished with the Space 
Telescope Imaging Spectrograph (STIS) and Far Ultravi- 
olet Spectroscopic Explorer (FUSE) spectrographs, which 
are sensitivity limited to the brightest QSOs and hence 
lowest rcdshifts. At these redshifts, few absorption lines 
are available and much of the determination of the phys- 
ical conditions in the systems is forced to be based on 
H I and O VI alone, or based on column density ratios 
of different ions such as C III and Si III. As such, non- 
solar abundance ratios of the different elements can be 
mistaken for the effects of ionization. 

In order to overcome some of these issues and to better 
constrain the physical conditions in the IGM and CGM, 
we have undertaken a blind survey of z ~ 1 QSOs with 
the COS spectrograph (Froning & Green 2009; Green et 
al. 2011). At moderate redshifts, we can access many 
of the numerous absorption lines in the FUV at A < 
912 A. Multiple transitions of other ionization states 
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AA 770,780 doublet to be observable. Other absorp- 
tion systems with multiple lines have been identified at 
z abs =0. 2552,0. 5778,0. 5910. 

We have used profile fitting to determine column densi- 
ties in the absorption lines detected in the PG1 148+549 
spectrum. The synthetic Voigt profiles were convolved 
with the COS line spread functions (LSF) 7 as deter- 
mined by the COS science team at the nearest tabulated 
wavelength (Ghavamian et al. 2009). Whenever possi- 
ble, multiple transitions were fit simultaneously (e.g. Ly- 
man /3, Lyman 7, Lyman S) to determine the best fitting 
column densities. Upper limits on column densities for 
non-detected species were determined in the manner de- 
scribed in Meiring et al. (2011). In short, synthetic Voigt 
profiles are created and the column density is increased 
until a 3cr limit is indicated given the S/N of the data. 

The COS wavelength solution is now known to have 
inaccuracies on one to two pixel levels ( i.e. up to ~ 30 
km s -1 ), especially at the wavelength extremes of the 
instrument (Savage et al. 2011). Shifts in velocity of this 
order have been applied to the absorption lines when 
necessary. The solar system abundances used throughout 
this paper have been adopted from Lodders (2003). 

3.1. The z a i, s =0.6838 system: 

Velocity plots for the lines are shown in Figure 1, with 
the best fit profiles overlayed in red. We have used two 
components to fit the absorption lines in this system. 
The Lyman a A 1215 line falls outside of the COS band- 
pass, however we detect transitions from Lyman f3 A 1025 
to Lyman e A 937. Further transitions in the Lyman se- 
ries are too weak to detect. 

Multiple metal lines are also detected. Lines of O IV 
A 787, O VI AA 1031, 1037, N IV A 765, and C III A 
977 are all detected at high significance. The O III A 
832 line is blended with the Milky Way Si IV A 1402 
line. By fitting the corresponding Si IV A 1393 line, we 
have removed this contamination from the O III A 832 
profile. A strong Lyman a line at z a b s =0.1529 is also 
observed at —175 km s _1 . This was also included in 
the fit to the O III line as some optical depth may be 
contaminating the profile due to the extended wings of 
the COS LSF. The fit to the O III A 832 line in Figure 
1 shows the contribution from Milky Way Si IV A 1402 
with the thick orange line. 

Comparisons of the apparent column density N a (v) 
profiles using the apparent optical depth method of Sav- 
age & Sembach (1991) are shown in Figures 2 and 3. The 
O IV and O VI profiles show excellent agreement in the 
component structure as can be seen in Figure 2. The 
apparent optical depth of the O IV A 787 line has been 
scaled by a factor of 0.75 for the sake of comparison. Col- 
umn densities via the apparent column density method 
and equivalent widths are given in Table 3, while the 
adopted column densities determined via profile fitting 
are given in Table 2. 

The Ne VIII A 780 line is blended with a C III A 977 
line at z o & s =0.5778. As this is the only C III line avail- 
able, it is impossible to remove the contamination from 
the Ne VIII A 780 line. The Ne VIII A 770 line is not 
blended with any other features, and we determine a col- 

7 Tabulated LSFs for COS are available at 
http : //www. stsci . edu/hst/cos/perf ormance/spectral_resolution/ 



umn density of log N(Ne VIII)=14.14±0.04 via profile 
fitting. The Ne VIII A 770 line is well aligned with the 
O VI A 1031 line, as can be seen in Figure 3, increasing 
our confidence in the veracity of the Ne VIII detection. 
This gives the ratio of N(Ne VIII)/N(0 VI)~0.42. 

3.2. The z a b s =0. 7015 system: 

No Lyman series absorption lines are detected in this 
system. We place an upper limit on N(H I) based on 
the non-detection of the Lyman j3 A 1025 line of log N(H 
I)<13.8. Velocity plots of the absorption lines from this 
system are shown in Figure 4, with the best fit profiles 
overlayed in red. 

Multiple metal lines are detected in this system. Lines 
of O III A 832, O IV A 787, O VI AA 1031,1037, C III 
A 977, and Ne VIII AA 770,780 are all detected. The 
Ne VIII AA 770,780 lines have a doublet ratio of 2:1, 
assuming the lines are weak and lie on the linear portion 
of the curve of growth the equivalent widths of the lines 
should have the same ratio. We find Wo(770)=28±5 mA 
and W o (780)=18±6 mA including errors from continuum 
placement and photon noise. The observed ratio of 1.6:1 
is close to the expected value, and there appears to be 
some contamination in the A 780 line in the red side of 
the profile which could account for the increased optical 
depth. The profiles do match well in terms of alignment 
in velocity space, as can be seen in Figure 7. 

Two components were used to fit the profiles at v = +3 
and v = +63 km s _1 relative to the systemic redshift 
of the system. The strongest component in each line is 
centered at v = +3 km s _1 . Apparent optical depth 
profiles of the O IV, O VI, and Ne VIII lines are shown 
in Figures 6 and 7. The O IV profile has been scaled by 
a factor of 1.5 for comparison to the O VI lines in Figure 
6. 

N III is not detected in this system, the N III A 989 
line is blended with a Lyman a line at z a f, s =0.3862, and 
the N III 685 line, which is nearly two times stronger, is 
not observed as well. The N IV A 765 line, the only N IV 
line available at these redshifts, is unfortunately blended 
with the strong Galactic absorption line of O I A 1302. 

With a simultaneous fit of the O VI AA 1031,1037 lines 
using two components in each line, we determine log N(0 
VI)=14.37±0.05 combined in the two components. The 
Ne VIII AA 770,780 lines were also fit simultaneously 
giving log N(Ne VIII)=13.86±0.05. This gives a ratio 
of N(Nc VIII)/N(0 VI)~0.30. Column densities via the 
apparent column density method and equivalent widths 
are given in Table 5, while the adopted column densities 
determined via profile fitting are given in Table 4. 

3.3. The z abs =0.7248 system: 

No Lyman series absorption lines are detected in this 
system. We place an upper limit on N(H I) based on 
the non-detection of the Lyman j3 A 1025 line of log N(H 
I)<13.7. Figure 5 shows multiple lines of interest in the 
system. 

The O IV A 787 and O VI AA 1031, 1037 lines are 
detected in this system at high significance. A single 
component centered at v = —4 km s" 1 relative to the 
systemic redshift of 2^=0.7248 was used to fit the pro- 
files. The adopted column densities determined via pro- 
file fitting are given in Table 6. 
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Figure 1. Velocity plots of the absorption lines in the z a ;, s =0.6838 system in PG1148+549. Vertical dotted lines mark the positions of 
the components. The solid red line represents the best fit models determined from profile fitting. The O III A 832 line is partially blended 
with the Si IV A 1402 line from the Milky Way. The contribution from Si IV A 1402 is shown with the solid orange line. Magenta point 
denote regions where the lines are blended with unaffiliated features from other systems. The cyan line at the bottom of each panel denotes 
the 1(T error of the normalized flux. 

A large suite of ions is necessary to constrain photo or 
collisional ionization models. As most of the O VI ab- 
sorbers observed to date are at lower redshift, the number 
of species available is limited to H I, O VI, Si III, C III 
and/or C IV in most cases. In a limited number of cases, 
a wider array of ions has been presented (e.g. Chen & 
Prochaska (2000); Sembach et al. (2004); Lehner et al. 
(2009)). 

The carbon to oxygen ratio is known to vary as a func- 
tion of metallicity, with carbon being sub-solar at the 
-0.5 dex level for [0/H]< -0.5 (Akerman et al. 2004; 
Pettini et al. 2008). The N V AA 1238,1242 lines are also 
sometimes seen in the IGM, and have also been used to 
constrain ionization models. The [N/O] ratio is also well 
known to be metallicity dependent in H II regions in the 
Milky Way and nearby galaxies (Van Zee et al. 1998; 
Van Zee & Haynes 2006). As such, the observed ratios 
of carbon to oxygen and nitrogen to oxygen in ionized 



Two features are located at the expected positions of 
the Ne VIII AA 770,780 lines with S/N~30 per resolution 
element in the region giving EW(770)=26±8 mA and 
EW(780)=19±5 mA Profile fitting of the lines gives log 
N(Ne VIII)=13.81±0.07, and a ratio of N(Ne VIII)/N(0 
VI) =0.89 indicating a high temperature if the gas is in 
collisional ionization equilibrium. 

The N IV A 787 and C III A 977 lines are heavily 
blended with features from other systems, and could not 
be fit. Apparent column density profile plots of the O 
IV, O VI and Ne VIII profiles are shown in Figures 8 
and 9. Column densities via the apparent column den- 
sity method and equivalent widths are given in Table 7. 



4. IONIZATION MODELING 
4.1. Photoionization 
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Table 3 

Equivalent widths and AOD column densities for transitions 
observed in the z a j,s=0.6838 system. 
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Figure 2. Apparent column density profiles of the O IV and O 
VI lines from the z ab3 =0.6838 system. The O IV A 787 line has 
been scaled by a factor of 0.75 for illustration. 
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Figure 3. Apparent column density profiles of the O VI and Ne 
VIII lines from the z abs =0.6838 system. The O VI A 1031 line has 
been scaled by a factor of 0.4 for illustration. 



Table 2 

Adopted column densities from profile fits to the COS data for 
the z aba =0.6838 system. 



Table 5 

Equivalent widths and AOD column densities for transitions 
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Table 6 

Measured column densities from profile fits to the COS data for 
the z a (, s =0.7248 system. 



"Due to the blending in this line, the b values for the components 
were fixed to the values determined from the O IV A 787 line. 
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Figure 4. Same as Figure 1, but for the 2 a (, s =0.7015 system. 
Table 7 

Equivalent widths and AOD column densities for transitions 
observed in the z a f, a =0.7248 system. 



Transition 


Wo 


N a (AOD) 




mA 


cm 


H I 1025 


<30 


< 13.7 


O II 834 


<15 


< 13.5 


O III 832 


<16 


< 13.5 


O IV 787 


27±5 


13.70±0.06 


O VI 1031 


71±20 


13.84±0.10 


O VI 1037 


37±12 


13.87±0.15 


Ne VIII 770 


26±8 


13.70±0.12 


Ne VIII 780 


19±5 


13.87±0.10 



plasmas are dependent on the metallicity as well as the 
temperature or ionizing flux, and make it rather difficult 
to tightly constrain parameters in the models. 

The optimal solution is to have a number of adjacent 
ions from a single species, such as oxygen for which there 
are a number of lines from different ionization states in 
the FUV. Nitrogen and sulfur also have a number of res- 



onance lines from multiple ionization states in the FUV, 
although the abundances of nitrogen and sulfur are lower 
than that of oxygen by a factor of ~ 7 and ^32, respec- 
tively, making the detection of the lines more challenging. 

To investigate the ionization mechanism in these sys- 
tems we have created grids of photoionization models 
using cloudy version 10.00, last described in Ferland et 
al. (1998). We have used the Haardt & Madau extra- 
galactic background spectrum (Haardt & Madau 1996), 
with the modifications and updates described in (Haardt 
& Madau 2001) at z = 0.7 including contributions from 
both AGN and galaxies to illuminate the plane-parallel 
clouds with ionizing radiation. 

We have matched the observed O III to O IV ratio to 
the models in order to constrain the ionization parameter 
U, which is the ratio of the number density of ionizing 
photons to the particle density U = rij/nn- The photon 
energy required to produce 3+ of 54.9 eV is just at the 
He + ionization potential of 54.4 eV. The atmospheres 
of O and B type stars produce photons at these ener- 
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Figure 5. Same as Figure 1, but for the z f, s =0.7248 system. 




Radial Velocity km s 1 Radial Velocity (km s ') 



Figure 6. Apparent column density profiles of the O IV and O 
VI lines from the z a ^ 3 = 0.7015 system. The O IV A 787 line has 
been scaled by a factor of 1.5 for illustration. 



Figure 7. Apparent column density profiles of the O VI and Ne 
VIII lines from the z ai)s =0.7015 system. The O VI A 1031 line has 
been scaled by a factor of 0.3 for illustration. 
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Figure 8. Apparent column density profiles of the O IV and O 
VI lines from the z aba =0.7248 system. The O IV A 787 line has 
been scaled by a factor of 1.5 for illustration. 
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Figure 9. Apparent column density profiles of the O VI and Ne 
VIII lines from the z ai)s =0.7248 system. The O VI A 1031 line has 
been scaled by a factor of 0.75 for illustration. 
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Figure 10. The extragalactic background spectra from Haardt 
& Madau (2001) in red, along with the modifications including an 
obscured AGN. Ionization potentials for several atoms included in 
the analysis of these absorbers are overplotted. 



gies, wheres at higher energies the photon flux steeply 
declines due to absorption from He + . Several lines of O 

II are covered in the data, although none are detected at 
adequate significance in any system and only weak upper 
limits can be placed on the columns density which do not 
constrain the models. 

In Figure 11 we show the results of the CLOUDY pho- 
toionization models. For both the z a h s =0.6838 and 
z Q &s=0.7015 systems the ionization parameter that best 
fits the O III and O IV column densities is log U=— 1.5. 
The O VI and Ne VIII column densities are severely un- 
derestimated at the ionization parameter best fitting the 
O III to O IV ratio. There is no solution that fits the O 
III, O IV, and O VI column densities simultaneously. O 

III is underproduced in the models matching O IV and 
O VI, while O VI is underproduced when O III and O 

IV are matched. 

As the O III A 832 line was not detected in the 
z Q &s=0.7248 system, we cannot place a unique value on 
the ionization parameter. We can however constrain the 
ionization parameter to a limited range. The observed 
ratio of O IV/O III>1.77 implies that log U> -1.75, 
while the observed ratio of O VI/O IV=1.29 bounds the 
ionization to be log U< —0.75. This allowed region is 
shaded in pink in the bottom row of Figure 1 1 . 

The close alignment between the low and high ions in 
these systems perhaps indicates a single-phase gas. As 
the Haardt & Madau (2001) background does not pro- 
duce the column density ratios observed in these systems, 
is it possible to modify the ionizing background spectrum 
to produce the amount of high ions observed in these sys- 
tems? 

The UV background is particularly uncertain in the 
soft X-Ray regime (100-500 eV), where these photons are 
easily absorbed by Galactic gas. In order to investigate 
whether modifications to the Haardt & Madau ionizing 
background could reproduce the observed column density 
ratios, we have also created similar grids of models with 
various spectral energy distributions (SEDs) in order to 
investigate whether alterations to the ionizing spectrum 
can produce the observed column densities. 

First, we have used the X-Ray emission from a so- 
lar metallicity 0.25 keV plasma, simulating the emission 
from the hot X-Ray emitting gas around nearby star 
forming galaxies which have been measured to have sim- 
ilar temperatures (Li et al. 2008). Such optically thin 
plasma does emit in the soft X-Rays, and is not included 
in the determination of the UV background of (Haardt & 
Madau 2001). This spectral energy distribution is dom- 
inated by emission lines, with some contribution from 
thermal continuum emission as well. 

The typical luminosity of an X-Ray emitting halo is 
~ 10 40 erg s" 1 (Li et al. 2008). We model the clouds as 
plane parallel slabs outside of this X-Ray emitting halo. 
We find that O VI and Ne VIII are not produced at the 
ratios observed in these systems until the flux reaches 
~ 2.5 x 10 -3 erg s _1 cm~ 2 , implying a distance from the 
source of ~ 0.2 kpc, assuming that the emission comes 
from a point source. In a real X-Ray emitting halo, the 
maximum flux is obviously lower as the emission arises 
in an extended volume. Nonetheless, such a small im- 
pact parameter is unphysical hence we rule this model 
out as a viable method to produce these high ions via 
photoionization. 
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We have also attempted to modify the Haardt & 
Madau (2001) background spectrum with the inclusion 
of an obscured AGN. With a total column density of 
~ 10 19 cm -2 , the spectrum of an AGN will have a cutoff 
at ~ 0.1 keV due to the photoelectric absorption of the 
neutral gas. This is just below the ionization potential 
of 4+ , and the low energy cutoff is introduced here was 
chosen specifically to generate a boosted flux in this en- 
ergy range. We have used a power law spectrum with 
the fiducial value of the spectral index a = 2.0. Includ- 
ing such a source will allow for extra photon flux in soft 
X-Rays which are capable of creating O VI and Ne VIII 
while keeping the photon flux low at longer wavelengths. 

Three models were created, with varying amounts of 
contribution to the soft X-Ray regime by increasing the 
intensity of the obscured AGN relative to the Haardt & 
Madau (2001) background. At a distance of 100, 500, and 
1000 kpc from an AGN with a typical luminosity of 10 45 
erg s" 1 , the AGN continuum will dominate the Haardt 
& Madau (2001) background by a factor of -1000, 40, 
and 10 respectively. These incident spectra are shown in 
Figure 10. 

Again, even with the increased flux at the ioniza- 
tion potentials for the high ions, the column density ra- 
tios predicted in the photoionization simulations do not 
match the observed values. As such, we are forced to 
assume that the gas has multiple phases. 

Could the high ions exist in a low density photoionized 
state? Such two phase photoionization models could po- 
tentially create a solution without a need for introducing 
a hot phase that is collisionally ionized. Problems arise 
however when investigating the sizes of the clouds in this 
scenario. The observed ratio of N(Ne VIII) /N(0 VI) is 
reached when log U ~ —0.25 for these systems, which 
in turn implies densities of — 1 x 10~ 5 and sizes of > 1 
Mpc. Such large structures would be affected by Hubble 
broadening of > 100 km s _1 , which is not supported by 
the best fitting Doppler parameters. Similar arguments 
have also been made for other systems with Ne VIII de- 
tections, where the implied cloud sizes in the case of pho- 
toionization production of Ne VIII lead to unrealistically 
large cloudy sizes (Narayanan et al. 2011). 

As single-phase, and even two phase photoionization 
models fail to account for both the low and high ions 
simultaneously, we are forced to conclude that the lower 
ions (O III, C III, O IV, etc.) reside in a photoionized 
core, and that the high ions (O VI and Ne VIII) reside 
in a collisionally ionized phase. Such necessity of a colli- 
sionally ionized phase where the high ions arise has also 
been demonstrated for other Ne VIII absorbers (Tripp 
et al. 2011; Narayanan et al. 2009, 2011; Savage et al. 
2005). 

4.2. Collisional Ionization 

As the photoionization models fail to produce the O 
VI and Ne VIII column densities observed, we include a 
hot collisionally ionized phased as well. We assume that 
gas is in collisional ionization equilibrium (CIE), and that 
the hot phase has the same metallicity as the photoionzed 
material. As the profiles of O IV, OVI, and Ne VIII are 
all well aligned in velocity space , the assumption that 
the hot gas has the same metallicity as the photoionzed 
gas seems reasonable. 

The steeply increasing Ne 7+ and steeply declining 5+ 



ionization ratios in CIE as shown in the right side of 
Figure 11 allow for a precise temperature of the hot gas 
to be determined for the systems investigated here. With 
the temperature determined from the ratio of Ne VIII to 
O VI, the total hydrogen column density is determined by 
scaling the N^* to match the observed column densities 
of Ne VIII and O VI. We find that log T-5.7 for these 
systems, similar to other absorption line systems with 
detected Ne VIII (Savage et al. 2005; Narayanan et al. 
2009, 2011; Tripp et al. 2011). 

We have also investigated the non-equilibrium colli- 
sional ionization models of Gnat & Sternberg (2007). 
Again, no single-phase model can account for the ob- 
served ratios of O III, O IV, O VI, and Ne VIII. Tem- 
peratures of log T~5.7 are also predicted in the non- 
equilibrium models based on the Ne VIII/O VI ratio. 

The Doppler parameters of the lines also provide a con- 
straint on the gas temperatures: b = y^2kT/m assuming 
the line profile is dominated by broadening from ther- 
mal motions in the cloud. This is effectively an upper 
limit of course, as turbulent motions can also contribute 
to the width of the profile. For log T^5.7 K as indi- 
cated by the Ne VIII/O VI ratios in these systems, the 
minimum Doppler parameter is bo vi > 22.7 km s~l or 
frjve vin > 20.6 to support the collisional ionization sce- 
nario. The Doppler parameters for all of the Ne VIII line 
in these systems are in agreement with the above con- 
straint. The best fit Doppler parameters for the O VI AA 
1031, 1037 lines in the z a b s =0.6838 and z Q f, s =0.7105 sys- 
tems are also fully consistent with collisional ionization. 
The z a bs =0.7248 system (which suffers from particularly 
noisy profiles) is slightly more discrepant, but still con- 
sistent within 2cr of the best fit Doppler parameter. 

5. RESULTS 
5.1. Physical Conditions 

As the modified background spectra described above 
do not reproduce the observed column densities, we pro- 
ceed under the assumption that the clouds are illumi- 
nated with the background radiation field of Haardt & 
Madau (2001). 

The extragalactic background spectrum of Haardt & 
Madau (2001) has an ionizing photon flux (hv > 13.6 
cV) of 2.2xl0 5 photons s^ 1 cm" 2 at z ~ 0.7. The back- 
ground intensity does vary slightly between z = 0.68 and 
z = 0.72, however for the ease of the calculations we as- 
sume the above value for all the systems. We estimate the 
density and cloud sizes of the photoionized portion of the 
gas as the ionization parameter U relates the photon flux 
(given by the ionizing background spectrum of Haardt & 
Madau 2001) to the density of the gas. We find that 
n ~ 10~ 4 cm~ 3 for these systems, corresponding to over- 
densities of A ~ 80 at z = 0.7 (A = p/pb)- Variations 
in the intensity of the ionizing background would affect 
the densities determined here. The ionizing photon flux 
4> is of course unknown, and almost certainly varies from 
place to place in the Universe with enhancements from 
local sources. It has been argued however that the lo- 
cal contributions are negligible for systems with log N(H 
I)<17.0 (Schaye 2006). Results for the photoionized por- 
tion of the absorbers traced by species such as H I, O III, 
C III and O IV are given in Table 8. 

We can predict the properties of the thermally broad- 
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Figure 11. From top to bottom, photo (left side) and collisional (right side) ionization modeling for the systems is shown for the 
• 2 abs = 0.6838, z a ;, s =0.7015, and z a [, s =0.7248 systems. For the 2 a j, s =0.6838 and z a f, s =0.7015 systems, the measured column densities are 
shown in large points along the vertical blue line which denotes the best fitting ionization parameter or temperature. For the z a (, s =0.7248 
system, the shaded pink area marks the allowed values of log U based on the O III, O IV, and O VI column densities. 



ened Lyman a lines in the collisionally ionized compo- 
nent in these systems implied by the CIE models de- 
scribed in §4. The Ne VIII/ O VI ratio gives a unique 
temperature for the hot phase of the gas bearing O VI 
and Ne VIII. The Doppler parameters for H I are pre- 
dicted to be 89-93 km s _1 at temperatures of log T~5.7. 
The thermally broadened Lyman a lines are likely to be 
difficult to detect for these systems for the following rea- 
sons. The optical depth at line center can be determined 
given the Doppler parameter and atomic data as 



r = 1.497 x 10 



_ 15 jV(cm- 2 )/A (A) 
b(km s _1 ) 



(1) 



For the z a b s ~0.6838 system, the thermally broadened Ly- 
man a line may be potentially observable, with a optical 
depth at line center of r = 0.38. This however is slightly 



misleading, as the Ne VIII and H I are well aligned in 
this system, so the broad feature will be superimposed 
on the much stronger Lyman a line from the photoion- 
ized gas. For the z a f, s =0.7015 and £ o f, s =0.7248 systems, 
the optical depth at line center is only r < 0.07; high 
S/N spectra may be able to able to detect such a fea- 
ture but small undulations in the QSO continuum make 
broad and weak features hard to detect. Parameters for 
the collisionally ionized phase of the absorbers are given 
in Tabic 9. 

5.2. Redshift Density and Baryonic Content 

The COS spectrum of PG1 148+549 enables detection 
of the rcdshifted Ne VIII AA 770,780 lines above z = 0.49, 
and we exclude a 3000 km s" 1 region at the emission red- 
shift of the QSO to not include any possible absorbers 
associated with the QSO itself (Tripp et al. 2008). As 
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the number of absorption lines in this spectrum is high, 
one must take into account the effects of line blanketing 
on the determination of the total absorption path length 
given in the equation above. Following the method out- 
lined in Tripp et al. (2008) to account for such effects and 
variable S/N, we find AX = 0.81 and Az = 0.41 for a 10 
pixel wide line and a minimum detection threshold of 20 
mA, corresponding to the smallest observed equivalent 
width for the weaker line of the Ne VIII doublet in these 
systems. 

We estimate the redshift density of Ne VIII absorbers 
to be dAf /dz~ 7±3 for W > 30toA given the three ab- 
sorbers detected in this spectrum over a total redshift 
path of Az — 0.41. We assume Poisson statistics in 
the determination of our confidence intervals, with val- 
ues taken from (Gehrels 1986). This estimate will be 
revised in future work when all of the spectra from this 
program will be analzed and presented. We note that 
the redshift density of O VI absorbers is only slightly 
greater than the estimate above for Ne VIII absorbers, 
with dAf OVI /dz~ 5 — 9 for O VI systems with equivalent 
widths similar to the systems here of EW(1031)> 80 mA 
(Tripp et al. 2008). 

We can make an estimate of the total comoving density 
traced by the O VI and Ne VIII absorbers from these 
observations. The comoving gas density is related to the 
column densities in the absorbers as 



a, 



Hp nm H J2" N i 
c p c AX 



(2) 



where AX is the absorption path length in a flat cosmol- 
ogy defined as 



AX = 



(l + zf 



s/ajl + zf + n A 



zdz 



(3) 



(Wolfe, Gawiser & Prochaska 2005). 

Using the N(Ht () values determined in Section 4 from 
matching the observed O VI and Ne VIII column den- 
sities in CIE, we find that the density of gas as traced 
by these absorbers is ft gas = 0.002, or roughly 4% of the 
baryon density determined by observations of the cosmic 
microwave background (Spergel et al. 2007). This is a 
crude estimate, and will be updated when the full sam- 
ple of Ne VIII absorbers in the program are analyzed. 
As ~ 50% of the baryons at low z are unaccounted for 
(Fukugita et al. 1998), there is still a sizable fraction of 
baryons missing. 



5.3. Sizes and Stability 
From the definition of the redshift density, 

dAf c 

= nirK — — 

dz H(z) 



(1 + ^) 2 , 



(4) 



the number density of absorbers n is directly related to 
the redshift density dAf /dz and cross sectional area of 
the absorbing cloud. Under the assumption that these 
systems are connected to the CGM in the halos of galax- 
ies, which seems likely as O VI systems have been clearly 
linked to the CGM of galaxies at lower redshifts (Tum- 
linson et al. (2011b); Prochaska et al. (2011), but see also 
Tripp et al. (2006)), we can estimate the typical sizes of 
the collisionally ionized clouds bearing O VI and Ne VIII 



from the observed redhsift density dAf /dz and number 
density of galaxies 4> at z ~ 0.7. 

The galaxy luminosity function is well constrained at 
the redshift of the systems in this work (Ilbert et al. 
2005). The survey of Ilbert et al. (2005) is reasonably 
complete down to 0.1L* at z = 0.7, as such, we take 
this as the lower bound in the integration over the lu- 
minosity function. The inclusion of lower luminosity 
galaxies would of course increase the overall number den- 
sity, however deep spectroscopic galaxy surveys of fields 
with known O VI absorbers indicate that dwarf galaxies 
(L< 0.L*) are unlikely to have a high covering fraction of 
O VI in their halos (e.g. Prochaska et al. 2011; Wakker & 
Savage 2009), and are hence unlikely to have an impact 
on the calculation. 

Taking the i?-band Schechter function parameters of 
M* = -22.62, <j) star = 1.9 x 10~ 3 , and a = -1.41 from 
(Ilbert et al. 2005), we find that n(L > 0.1L*) = 2.4 x 
10~ 2 . This is turn yields characteristic sizes for the O 
VI+Ne VIII systems of ~ 70 — 150 kpc. Such sizes are 
quite similar to what has been determined for the sizes 
of the halos of lower redshift galaxies via absorption line 
statistics (Tumlinson et al. 2011b; Prochaska et al. 2011; 
Wakker & Savage 2009). 

For the photoionized portions of the clouds, wc can es- 
timate the stability to collapse as we have both a density 
and size from the photoionization modeling. The Jeans 
length gives the characteristic length scale for which a 
cloud can be stable against gravitational fragmentation 
and collapse. Given a gas fraction / s , density nni, and 
temperate T, the Jeans length is given by (Schaye 2001) 



Lj ~ 6.5 x 10 



4 pc(- 



1/2 / t 



1/2 



(5) 

The CLOUDY simulations described in §4 determined 
temperatures of T ~ 2 x 10 4 K for the photoionized gas, 
combined with the densities given above of rag ~ 2 x 10 -4 
gives Jean's lengths of Lj ~ 65 kpc assuming that the 
clouds are gas dominated. All of the systems in this 
work are constrained to have sizes well below Lj, imply- 
ing that the ultimate fate of these systems is to evapo- 
rate into the surrounding medium (unless they are con- 
fined by some other mechanism). Similar such metal rich, 
compact clouds have also been observed at high redshift 
(Schaye et al. 2007). 

6. AFFILIATED GALAXIES 

We have obtained deep multiband imaging of the 
PG1148+549 field with the Large Binocular Telescope 
(LBT) and Large Binocular Camera (LBC). The U-band 
image of the field centered on the QSO is shown in Fig- 
ure 12. The data have a limiting magnitude of mu < 27. 
Spectra of galaxies labelled 76_24 and 142_30 were ob- 
tained with the Keck Low Resolution Imaging Spectrom- 
eter (LPJS). The data were reduced and analysed as de- 
scribed in (Werk et al. 2011). Parameters for the galaxies 
labelled in Figure 12 are given in Table 10. 

Strong absorption from C III A 977, O VI AA 1031, 
1037, and multiple Lyman series lines are seen at 
z a b s =0.2525 associated with galaxy 76_24. More inter- 
esting for this work, the galaxy 142_30 is only Av ~ 
+30 km s -1 from the absorber at z a bs= 0.7248. This 
galaxy is bright, with L^L*, has a high metallicity 
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Table 8 

Conditions in the photoionized gas of the z 



0.7 absorbers. 



Redshift 


log N(H I) 


log N(Htot) 


log U 


[X/H] 


Density 


Size 




cm -2 


cm z 






cm~ 3 


kpc 


0.6838 


14.79±0.03 


18.80 


-1.5 


-0.5 


2xl0" 4 


10 


0.7015 


< 13.8 


<17.69 


-1.5 


> 0.0 


2xl0 -4 


< 1 


0.7248 


< 13.7 


< 18.46 


-1.75 to -0.75 


> -0.5 


(0.4-4.0) xl0~ 4 


< 23 



Table 9 

Parameters for the collisionally ionized gas of the z ~ 0.7 
absorbers. 



Redshift 


logT 


log N(H tot ) 


log N(H I)« 


b 6 A 

prea 


~0 




K 


_9 

cm 


_9 

cm 


km s 




0.6838 


5.68 


19.80 


13.65 


89 


0.38 


0.7015 


5.69 


19.15 


12.90 


90 


0.07 


0.7248 


5.72 


18.83 


13.08 


93 


0.10 



"The predicted neutral 
temperature and N(Ht D t). 
'The predicted Dopplcr 
temperature. 
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Figure 12. A deep U band image of the PG1148+549 field cen- 
tered on the QSO. Galaxies for which we have obtained spectro- 
scopic redshifts are labeled according to their position angle and 
distance in arcseconds from the QSO. 

of [O/H]=+0.22±0.15, and a relatively high SFR of 
6.42±0.44 M Q yr _1 . The LRIS spectrum of this galaxy 
is shown in Figure 13. 

This high metallicity is consistent with the constraints 
on the metallicity of the absorption system of [0/H]> 
0.0. Similarly, the high star formation rate is consis- 
tent with the recent findings of Tumlinson ct al. (2011b) 
and Prochaska et al. (2011) that O VI is ubiquitously 
observed in the halos of star forming galaxies. Several 
other galaxies are also clearly seen in Figure 12. A more 
complete redshift survey of these galaxies is required to 
eliminate the possibility of other galaxies at the absorp- 
tion redshift as well. 
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Figure 13. Keck/LRIS spectrum of the galaxy located at z em = 
0.7248 at an impact parameter of 217 kpc. Prominent spectral 
features are labeled. The cyan line below the spectrum represents 
the la error in the flux. 

7. SUMMARY AND DISCUSSION 

In this work we have reported on the physical condi- 
tions and ionization in the Ne VIII + O VI absorbers in 
the COS spectrum of PG1 148+549. We summarize the 
results as follows: 

1. A single-phase model in which all of the ions are 
produced via photoionization or collisional ioniza- 
tion fails to reproduce the observed column den- 
sity ratios. Models in which the high ions are 
produced via photoionization predict unphysically 
large cloud sizes, with absorption lines that would 
be dominated by Hubble broadening. We conclude 
that the O VI and Ne VIII arise in collisionally ion- 
ized gas with T~ 10 5 7 K, while the low ions are 
produced via photoionization. 

2. The metallicities of the systems, as determined 
from the low ions and photoionization modeling, 
are all [0/H]> —0.5, higher than the typically as- 
sumed metallicity of the IGM of [X/H] 1.0. The 

z a bs=0-7015 and z a fc s =0.7248 systems indicate su- 
persolar metallicities. 

3. We have investigated the consequences of modify- 
ing the ionizing background of (Haardt & Madau 
2001). We have attempted to include an obscured 
AGN with a soft X-Ray cutoff at 100 eV, which 
would boost the flux near the ionization potentials 
required to produce O VI and Ne VIII. Even with 
these significant modifications to the extragalactic 
background of Haardt & Madau (2001), the mod- 
els fail to produce the observed O VI and Ne VIII 
column densities. 

4. We have determined that the gas as traced by O 
VI and Nc VIII has temperatures of T~ 10" K 
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via the VI to Ne VIII ratio in collisional ioniza- 
tion equilibrium. Non-equilibrium models of Gnat 
& Sternberg (2007) also give similar temperatures 
based on the O VI to Nc VIII ratio. 

5. Densities of nn ~ 10~ 4 cm" 3 , corresponding to 
overdensities of A ~ 80 at z = 0.7 arc deter- 
mined from the photoionization modeling. The 
photoionzcd portions of the clouds are compact, 
with sizes R< 10 kpc. 

6. The Jean's length for the photoionizcd portion of 
these systems are Lj ~ 65 kpc, much larger than 
the derived cloud sizes. The clouds are thus not 
gravitationally bound, and will evaporate while en- 
riching the IGM and CGM with metals. 

7. We determine the redshift density of Ne VIII ab- 
sorbers to be cW /dz~ 7+3 for Wo > 30mA, with a 
cosmic baryon density of fl gas = 0.002, or roughly 
4% of the baryon density determined by observa- 
tions of the cosmic microwave background. A large 
fraction of the baryons are thus still unaccounted 
for. 

8. Based on the observed redshift density and galaxy 
luminosity function at z ~ 0.7, we estimate the 
characteristic size of the collisionally ionized gas 
harboring the O VI + Nc VIII absorbers to be ~ 
70 - 150 kpc. 

9. Wc find a luminous star forming galaxy at an im- 
pact parameter of ~ 220 kpc from the z a bs =0.7248 
system. The galaxy also has a high metallicity, 
linking the absorpion line system to the hot ex- 
tended halo of the galaxy. 

It is interesting that the profiles of the O III, O IV, and 
O VI are well aligned (in the sense that the component 
structure is similar for all the lines), but a single-phase 
model consisting of only photoionization or collisional 
ionization fails to reproduce the observed column densi- 
ties. Based on the ionization modeling we are forced to 
conclude that the highly ionized species of O VI and Nc 
VIII are not created in by pure photoionization in these 
systems. 

Turbulent mixing layers are another mechanism in 
which to create highly ionized ions in the halos of galax- 
ies. In this scenario, cool, dense material is mixed with 
hot (~ 10 6 K) gas as the materials slide past each other 
(Begelman & Fabian 1990; Slavin & Cox 1993; Kwak & 
Shclton 2010). Kelvin- Hclmholtz instabilities in the gas 
create a mixing of the phases and the mixed gas produces 
elevated amounts of high ions than collisional ionization 
equilibrium (Kwak & Shelton 2010; Kwak et al. 2011). 

Recent COS observations indicate a near unity filling 
factor for O VI within 150 kpc of star forming galaxies 
(Tumlinson et al. 2011b). If we assume a typical cloud 
size of 10 kpc similar to what we have determined here, 
^200 individual clouds would be necessary along a given 
line of sight to ensure that the covering fraction is near 
unity. The small cloud sizes determined here for the pho- 
toionzcd gas as traced by species such as H I, O III, and O 
IV would require many clouds in order to produce a high 
filling factor such as is seen in (Tumlinson et al. 2011b). 
If the O VI arises in a thin skin on the surface of the 
lower-ionization clouds, then the size of the O VI cloud 
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ID p SFR W r [O/H] 
kpc Mq yr — 1 

76.24 0.25250 95 1.53±0.04 -19.86 ±0.04 -0.11 ±0.15 

142.30 0.72492 217 6.42±0.44 -21.15 ±0.13 +0.22 ±0.15 



would be effectively similar to that of the low-ionization 
phase. Such a scenario would also likely lead to much 
more complex profiles than what is typically seen (some 
systems do however show remarkably complex profiles in 
VI however Tumlinson et al. (2011a)). 

The turbulent mixing layer scenario, while it does pro- 
duce elevated levels of high ions, would still require many 
clouds in order to account for the large filling factor of 
O VI seen in Tumlinson et al. (2011b) as the high ions 
are produced on the skins of the clouds. One possible 
explanation is that the O VI is resulting from the evap- 
oration of the outer parts of the clouds expanding into 
the halo of the galaxy, leaving behind a smaller and less 
ionized cloud that is traced by H I, Mg II, C II, etc. 
It is also unclear if the turbulent mixing layers would 
be effective at large impact parameters to the galaxies, 
where the densities are low and mixing would therefore 
less effective. 

Much effort has been made recently into simulations 
of structure formation including the enrichment of the 
IGM via galaxy feedback with hydrodynamical simula- 
tions (Tepper-Garcia et al. 2011; Dave et al. 2001; Op- 
penheimer & Dave 2006, 2009; Oppcnheimer et al. 2011). 
The origins of the O VI absorbers is of much debate. 
Some simulation indicate that the O VI absorbers are 
mainly collisionally ionized (e.g. Cen & Chisari 2011), 
while others point to photoionization as the main mech- 
anism producing O VI (e.g. Oppenheimer & Dave 2009; 
Oppcnheimer ct al. 2011). Some simulations also indi- 
cate that the stronger systems with EW> 100 mA are 
more likely to be collisionally ionized (Cen et al. 2001; 
Fang & Bryan 2001; Oppenheimer & Dave 2009; Op- 
penheimer et al. 2011). The observations presented here 
support a collisional ionization mechanism for the O VI 
and Ne VIII, in general agreement with the hydrody- 
namic simulations as the equivalent widths of the O VI 
lines are among the stronger systems that arise in the 
simulations. 

The metallicities of O VI absorbers are similarly dis- 
parate in the simulations. The simulations of Tepper- 
Garcia et al. (2011) and Oppenheimer & Dave (2009) 
indicate that most of the O VI absorbers reside in metal 
enriched regions with —1.0 <[Z/Zq]<0.0, whereas in Cen 
& Chisari (2011) the metallicities of the O VI absorbers 
is much lower, with —2.5 <[Z/Zq]< —0.5. The metal- 
licities determined here, as well as other systems in for 
which there arc multiple ionization states to adequately 
constrain the ionization models similarly indicate metal- 
licities of [0/H]> -0.5 (Howk et al. 2009; Narayanan ct 
al. 2011; Savage et al. 2005; Tripp ct al. 2011). 

Based on a large sample of QSOs, Tripp et al. (2008) 
find a redshift density of dAA /d2~ 5 — 9 for systems 
with Wo > 100 mA. The similar value determined for the 
Ne VIII systems of &Af NeVln /dz~ 7 perhaps indicates 
that many O VI systems could harbor Nc VIII lines. 
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Ne VIII may only be detectable in the strongest O VI 
absorbers as the neon abundance is lower than that of 
oxygen by a factor of ~7 (Lodders 2003). 

As was noted in the previous sections, the exact shape 
of the ionizing background is unknown and a limiting fac- 
tor when determining the mctallicitics and physical con- 
ditions in QSO absorption systems. With multiple ions 
of the same element, or perhaps multiple ions of several 
elements which have little intrinsic scatter in their rela- 
tive abundances, it may be possible to measure the shape 
of the ionizing background radiation directly based on 
the relative proportions of column densities in different 
ionization stages. Oxygen and sulfur could be possibly 
be utilized simultaneously as both are a-elements with 
multiple ionization stages present in the FUV. We plan 
to investigate this possibility further in future work. 
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